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ABSTRACT 



O Context. Asymptotic Giant Branch (AGB) stars are major contributors to both the chemical enrichment of the interstellar medium 
and the integrated light of galaxies. Despite its importance, the AGB is one of the least understood phases of stellar evolution. The 
main difficulties associated with detailed modelling of the AGB are related to the mass-loss process and the 3'^'' dredge-up efficiency 
Aims. To provide direct measures of mass-loss rates and luminosities for a complete sample of AGB stars in the Large Magellanic 
Cloud, disentangling the C- and O-rich stellar populations. 

Methods. Dust radiative transfer models are presented for all 374 AGB stars candidates in one of the fields observed by the new 
^ . VISTA survey of the Magellanic Clouds (VMC). Mass-loss rates, luminosities and a classification of C- and O-rich stars are derived 

5— ( ' by fitting the models to the spectral energy distribution (SED) obtained by combining VMC data with existing optical, near-, and 

' mid-infrared photometry. 

, Results. The classification technique is reliable at a level of - at worst - 75% and significantly better for the reddest dusty stars. We 

classified none of the stars with a relevant mass-loss rate as O-rich, and we can exclude the presence of more than one dusty O-rich 
star at a ~ 94% level. The bolometric luminosity function we obtained is fully consistent with most of the literature data on the LMC 
and with the prediction of theoretical models, with a peak of the C-star distribution at Mboi - -4.8 mag and no stars brighter than the 
^ ' classical AGB tip, at Mboi = -7.1 mag. 

Conclusions. This exploratory study shows that our method provides reliable mass-loss rates, luminosities and chemical classifi- 
cations for all AGB stars. These results offer already important constraints to AGB evolutionary models. Most of our conclusions, 
nI ' especially for the rarer dust-enshrouded extreme AGB stars, are however strongly limited by the relatively small area covered by our 

\l I study. Forthcoming VMC observations will easily remove this limitation. 

I Key words, stars: AGB and post-AGB - stars: mass loss - Magellanic Clouds 

1. Introduction noring crucial aspects of the AGB evolution. A more success- 
Zl ' rr., ^ ■ ^- , ,»^T,s • , , ful approach is based on synthetic codes, in which the stel- 
^ • The Asymptotic Giant Branch (AGB) is the last stage of ac- evolution is described by means of simplified relations de- 
^ , tive nuclear burnmg for low and intermediate mass (-0.8-8 M^) ^-^^^ ^^^^ complete stellar models, while convective dredge- 
d ■ stars. Although it i s a short-lived phase - less than a few Myr m ass loss are tuned by means of a few adjustable pa- 
(| GirardietalJ | 2010|) - AGB stars are major polluters of the in- ^ameters (iRenzini & VolilfTosl iGroenewegen & de Jo^fTool 



terstellar me dium and rnajor contribu tors to the integrated Ught ^fari^|2002HCordier et al.l 120071) . This approach provided the 



of galaxies dRenzini & Buzzonj 19861) In spite of its relevance, ^^^^ ^^^^^^^^ ^^^jl^^ isochrones suitable for reproduction of 

the AGB IS one of the least understood phases of stellar eyolu- ^^^.^ ^^^^^^^^ properties o f AGB stars, a nd in particular the 

tion, and a major source of disagreement between the results of ^^- ^^ ^^^^^ ^^^^^ (iMarigo et al.ll2008h . Despite the 

different population synthesis rnodels. This is mainly due to the ^^^^^^^ ^^^^^ ^^^^1^^ ^^^^ ^^jU ^^^^ parameters 

difficulties associated with modelling convective dredge-up and 3^ dredge-up and mass-loss efficiency) that need a detailed 

mass-loss processes. calibration, as shown by the disagreement between the observed 

Most of the available sets of evolutionary tracks and ^^^^^^ luminosity function (LF) o f C- and O-rich s tars in 

isochrones cover AGB ev olution in a very approximate way, ig- ^earbv gal axies and model pre dictions dGullieuszik et allM 

~in^ — : Z 7- J .r^o^ J iHeld et al.ll2010t iGirardi et al.ll2010l) . A more detailed calibra- 

Based on observations made with VISTA at ESO under program ID . ^ , f \ : ^ , , 

179 B-2003 models requires measures of mass-loss rates and a 
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reliable classification of C-rich and O-rich AGB stellar popula- 
tions in nearby stellar systems. 

Direct measures of mass-loss rates in AGB stars are lim- 
ited to samples of a hundred stars in the Mag ellanic Clouds 
dvan Loon et al.ll2b08t iGroenewegen et aDl2009l) . Estimates us- 
ing larger samples of stars are, to date, based on empirical re- 
lations betwe en mid-infrared (IR) colours and mass-loss rates 
(iMatsuura et al . 2QQ9). Recently, iSrinivasan et al.l (l2009l) used 



mid-IR photometry from the SAGE survey ( iMeixner et al.i2006l) 
to identify thousands of AGB star candidates in the LMC. Their 
spectral energy distributions (SED) were fitted with dust-free 
photospheric models for O- and C-rich stars. The mass-loss rates 
were then derived from the excess in the fluxes observed in the 
8 and 24 fim bands. This was the first step of an ongoing project 
aimed at directly measuring mass-loss rat es for all their cand i- 
dates using a grid of dust shell models (see lSargent et al.l[201 ih . 

To calibrate the 3'^'' dredge-up efficiency complete catalogues 
of O- and C-rich AGB stars are needed. The most straightfor- 
ward method to obtain a reliable classification uses a spectro- 
scopic classification, but this would be extremely expensive in 
terms of observing time. The most recent and impress ive effort 
in this direction is the SAGE-Spec d Woods et al.ll201 ll) spectro- 
scopic survey, which aims at a classification of mid-IR sources 
(mainly AGB stars) in the LMC. The project was awarded 224.6 
hr of Spitzer Space Telescope observations, and provided a reli- 
able classification of 197 point sources. This is clearly a minor 
fraction of the total LMC AGB population, whic h consists of 
some tens of thousands of stars dCioni et alJl2006h . A more ef- 
ficient, although in principle less reliable way to achieve a clas- 
sification employs photometric data. Different approaches have 



near-IR broad-band photometry (e. 


2.. Battinelh & Demersl2005l 


iGullieuszik et al.ll2008l;lHeld et al.l 


2010h. 



The Magellanic Clouds (MC) are ideal targets to study the 
details of AGB evolution, since they host a huge population 
of AGB stars covering the whole age (i.e. initial mass) range 
and they are nearby galaxies, easily accessible for observations 
of resolved stellar populations. The immediate drawback of the 
proximity is the large sky area covered by the MCs. Only re- 
cently have some observational campaigns provided coverage of 
a significant fraction of both the LMC and the SMC, in particu- 
lar: the optical Magellanic Clouds Photometric Survey (MCPS, 
IZaritskv et all2002ll2004 : t he all-sk y 2MASS dSkrutskie et al.l 
2006) and the MCs IRSF dKatp_et al. 2 (307) s urveys in the 
near-IR; the SAGE -LMC (Meixner etal] l2006l) . SAGE-SMC 
dGordon et al.ll201lh . S^MC tBolatto et aljIgOOTh Spiter surveys. 



and the AKARI/IRC LMC survey dita et al.ll2008h in the mid-IR 

The VISTA Magellanic Cloud survey (VMC, ICioni et all 

1201 ih is one of the six ESO public survey projects 
that are carried out with th e new ESO telescope VISTA 
dEmerson & Sutherland! I20T0I) . It aims at imaging about 180 
square degrees in the Magellanic system (LMC, SMC, the 
Bridge and the Stream) in the YJKs wavebands, reaching a depth 
of Ks ~ 20.5 mag, nearly 5 mag deeper than 2MASS and ~ 3 
mag deeper than IRSF. Observations started in November 2009 
and will take about five years to reach completion. The first im- 
ages and catalogues were released to the community in 2011 
and new data will be released at regular intervals according to 
the ESO policy for public surveys. The VMC data will provide, 
e.g., a detailed history of star formation across the Magellanic 
system and a measurement of its 3D geometry. The VMC is 
described in more detail in ICioni et al . (2011); a series of pa- 
pers will prese n t the scientific results of the survey. To date, 
iMiszalski et al.l d201 ih presented a multi-wavelength study of 



LMC planetary nebulae and lRubele et al.l (1201 ih obtained a spa- 
tially resolved Star Formation History of the LMC fields already 
observed by the VMC. 

We aim at extending the direct measure of luminosity and 
mass-loss rate to the whole AGB stellar population in the 
Magellanic system by modelling the SEDs obtained by collect- 
ing all available photometry from optical to mid-IR wavelengths. 
VMC observations are particularly important, because their pho- 
tometric quality and spatial resolution are much higher than that 
of 2MASS. This yields a more accurate photometry, in partic- 
ular for the extremely red dust enshrouded AGB stars, around 
the detection limit of 2MASS in the bluest bands - J ~ 16.5 
mag and K, ~ 15.0 mag. In addition AGB stars are variables 
and the addition of VMC data -taken at a different epoch from 
e.g. 2MASS- will result in a more robust estimate of the mean 
luminosity and mass-loss rate of the star The SED analysis also 
provides a classification of the stellar chemical composition. To 
test its power to disentangle O- and C-rich AGB stellar popula- 
tions, we present a critical comparison between our results and 
spectroscopic classifications from the literature. 

In this paper we present our first results, based on data of the 
LMC region covered by the first VMC observations. The data 
and the AGB selection criteria are presented in Sect. |2j Section 
|3] describes our model; the results are presented in Sect. |4l in 
Sect. |5] we summarise our work. 

2. The data 

We used the vl.O VMC release, which is described in detail in 
ICioni et all d2()l ll). The data were p rocessed b y the VISTA Data 
Flow Systern femerson et al.ll2004l) pipeline dirwin et al. ll2004h 
and r etrieved from the VISTA Science Archive ( Hambl v et al.l 
120041) . The major issue for AGB stars in VMC photometry is 
saturation and non-linearity of the detectors for bright objects. 
The VISTA pipeline provides a correction for saturation effects 
which can recover reliable photometry for stars up to ~ 2 mag 
brighter than the saturation limit. This is however not enough, 
since some AGB stars are even brighter, so we had to reject VMC 
photometry for stars brighter than 10.5 mag in the Y and / bands 
and 10.0 mag in the band (see lCioni et al ] |20TT1) . 

Additional photometry was obtained by collecting data from 
LMC photometric surveys: 



- optical: UBVI data from the MCPS (IZaritsky et al.ll2004 ; 

- near-IR: JHK , data from flie all-sky 2MASS 

dSkrutskie et al.l l2006l) . and the extended mission 6x 
long-exposure release; stars with very poor photometry - 
quality flag "E" - in any of the JHKs bands were discarded; 

- mid-IR: IRAC 3.6, 4.5, 5.8, 8.0 and MIPS 24, 70, 16 urn 
data from the SAGE catalogue dMeixner et al. I I2OO6I) : we 
used the SAGELMCcataloglRAC and SAGEcatalogMIPS24 
catalogues containing both SAGE epochs 1 and 2 separately 
to avoid losing information on variability; 

- mid-IR 9 and 18 fim ph otometry froni the A KARI/IRC mid- 
infrared all-sky survey dishihara et alJl2010h . 

The SED of each star was obtained using all photometric 
data points. When a single magnitude for each photometric band 
was needed, e.g., to plot colour-magnitude diagrams (CMDs), 
the data were merged. For mid-IR photometry we used the orig- 
inal SAGE epoch-merged data, while VMC and 2MASS / and 
Ks were averaged after applying the colour equations to trans- 
form J and Ks 2MASS p hotometry to the VISTA photometric 
system dCioni et al .11201 11) . 
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For our project the use of optical and mid-IR photometry is 
fundamental, so we will limit our analysis to the areas of the 
MCs observed by all of the VMC, the MCPS and Spitzer To 
date, there are only two VMC fields in common wit h the optical 
and m id-IR surveys, namely fields 6_6 and 8_3 (see ICioni et al.l 
The first is centred on the 30 Doradus star-forming re- 
gion, which is affected by severe crowding and differential red- 
dening problems. Since in this paper we want to define our 
methodology and test its reliability, it is clear that data obtained 
from observations of field 6_6 are not a suitable benchmark. We 
consequently base the analysis presented in this paper on data in 
VMC field 8_3, and in particular in the 1.42 deg^ region delin- 
eated by 75°50 < a(J2000) < 7T.5Q and -66?95 < (5(J2000) < 
-65° 19. 

The targets for our analysis were selected using near- and 
mid-IR photometry. We used 2MASS magnitude for selection 
rather than VMC ones, since some bright AGB stars have no 
reliable VMC photometry because of saturation. 

The only problem with this choice could be those AGB stars 
that are heavily obscured by thick circumstellar dusty discs with 
extremely red colours that are too faint in the near-IR for detec- 
tion by 2MASS. We verified this hypothesis and checked that 
no such sources exist in the current sample. The conditions we 
imposed are: 

- J - Ks > -0.07 5 X K, + 1.85, the relation defined by 
ICioni et al] ( |2006|) to minimise contamination by foreground 
Milky Way and LMC red supergiant stellar populations. 

- Ks brighter than 1 2.0 mag, which c orresponds to the tip of 
the RGB (TRGB) dCioni et al.ll2000h . This condition was ap- 
plied only to stars with J - Kg < 1.5 mag, since redder stars 
cannot be RGB stars, and extremely red dust enshrouded 
AGB stars can be fainter than the Ks-hand TRGB. 

- IRAC 4.5 fim magnitude brighter than 12.0 mag. This con- 
dition was set to avoid sources with faint mid-IR emission 
which cannot be AGB stars, but may be rather artifacts in 
the 2MASS catalogue. 

These criteria define a sample of 372 candidate AGB stars. We 
note that only 15 of them have 9 yL/m AKARI photometry and 
only 4 are detected in the 18 fim AKARI band. None are present 
in the AKARI/FIS All-Sky Survey Bright Source Catalogue. 
Finally, only one of the selected sources was detected in the 
MIPS 160 /im band and three at 70 jum. All of these objects 
turned out to have SEDs that are incompatible with AGB stars 
(see next section). 

The selection criteria defined above could in principle ex- 
clude th e most extremely red AGB stars, like those discov- 
ered by iGruendl et all (I2008h using mid-IR photometry and 
spectroscopy. Their SEDs peak in the mid-IR at wavelengths 
> IQfim and they are extremely f aint in bluer wave-bands. 
Some of the stars in iGruendl et al.l (l2008h sample have IRAC 
4.5 jim magnitudes fainter than 12.0 mag, and they could be 
likely fainter than the VMC detection limit. For this reason, we 
added an additional selection criterion, based on mid-IR pho- 
tometry, to include extremely red stars in our sample. U sing all 
the sta rs classified as candidate AGB stars by Gruend l & Chul 
(l2009l) . we defined the following selection rules: [4.5] -[8.0] > 1 
and [8.0] < 8.5 + 0.25 x ([4.5] - [8.0]), which correspond to the 
box displayed in the lower-right panel of Fig. [1] This criterion 
adds two objects to our sample. Anticipating the results that will 
be presented in the next sections, we can safely exclude that the 
bluest of the two objects is an AGB star. The classification of 
second one is quite unclear, but in the Appendix of this paper we 



provide some evidence that also this object could not be an AGB 
star 



3. The model 

We fitted the SEDs with a combination of photospheric mod- 
els, dust models and by solving the radiative transfer equa- 
tion. T he code used in this p aper is based on that presented 



by Groenewegen et al 



- (I2009I) . updated to use the DUSTY 

code dlvezic et alj 19991) to solve the radiative transfer equation 
(Groenewegen, in pr ep.). Another important update is the use of 
lAringer etal1(l2009h models for the stell ar atmosphere of C-rich 
stars. For O-rich stars MARCS models jGustafsson et al.ll2008l) 
were used. 

We a dopted only two di f ferent dust-grain compositions, 
following iGroenewegen et all (|2009|) . A 94% amorphous car- 
bon + 6% SiC mixture was assumed for C-rich stars, and 
"astronomical silicates" with the absorption coefficients of 
Volk & Kwqk (1988) for O-rich stars. Since we are considering 
on ly broadband photometry, we preferred to apply the results 
of IGroenewegen et alJ (l2009l) . based on mid-IR spectroscopy, 
rather than introducing the dust composition as a free parame- 
ter weakly constrained by the data used in our analysis. 

The observed SED is fitted using a minimisation technique 
to find the best-fitting values for the luminosity and optical depth 
T,i at 0.55 jum of the dusty shell under the assumption of a 
density profile. The minimisation is based on reduced (Xi)^ 
which was calculated by taking into account the photometric 
er rors. The mass-loss r ate M was computed using the relation 
of IGroenewegen et al.l (Il998h . applied to a geometrically thick 
shell, since we assumed an outer to inner radius ratio of 10"* 
(IGroenewegen et aDl2009l) : 



r.) = 5.405 X 10*^ 



rd R* Vexp Pd ' 



(1) 



where M is in Mq yr Vexp is the shell expansion velocity in 
km s \Ri, is the stellar radius in solar radii, is the inner dust 
radius in stellar radii, 2.1 is the absorption coefficient, a is the 
dust grain radius in cm, 'I' is the dust-to-gas mass ratio, and pd 
is the grain density in g cm""*. The values of Q^/a and Ri, were 
obtained from the dust and stellar models, respectively, while 
is an output of the dust radiative code. We adopted Vexp = 10 km 
s ' and *!* - 0.005. These are standard values for AGB stars in 
our Galaxy, which may not apply to the m ost metal-poor stars in 
the LMC (see, e.g.. lvan Loon et al.ll2008l) . 

The sample of stars was divided into a "Red" and a "Blue" 
group, using as a discriminant the condition J - > 1.5 mag. 
The number of stars in the red and blue groups is 97 and 274, 
respectively. The first group consists of the so-called red tail, the 
location in the CMD where most of the mass-losing C-stars are 
found. With this choice, in the "Red" group there are mostly low 
effective temperature dusty stars, whose SEDs are dominated by 
dust emission. In this case the choice of the atmospheric model 
is therefore not critical. On the other hand, bluer stars are ex- 
pected to have higher T^ff and lower mass-loss rates and their 
SEDs are determined largely by the stellar photosphere rather 
than the dust envelope, making the choice of photospheric model 
more critical. For the stars in this group we used a wider temper- 
ature range. Since the bluest stars are expected to have very low 
to negligible mass-loss rates, for all temperatures we considered 
the possibility of having dust-free models, fixing the mass-loss 
rate at zero. Only for the two lowest temperatures we left the 
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Fig. 1: Near- [upper panels) and mid-IR (lower panels) CMDs for stars in the VMC field 8_3. Stars classified as O-rich (upper left 
panel, blue dots) and C-rich (upper right panel, red dots) are shown as solid symbols. The five non-AGB sources are displayed 
as diamonds. In the upper panels we plotted curves of constant b olometric magnitudes (M hni between -6.0 and -4.0, with steps 
of 0.5 mag) derived using the A'j-band bolometric correction from lKerschbaum et al.l (|2010|) . In the lower-left panel all AGB stars 
are shown as bigger points, while in the lower-right panel bigger (green) points show all stars not classified as AGB candidates 
using the near-IR crite ria. The small grey cro sses and open circles show all the sources in the LMC classified as AGB stars and 
EROs, respectively, bv iGruendl & Chul (|2009|) . The filled circles represent the two additional stars selected purely on the basis of 
the mid-IR photometry using the area delineated by the solid line. 



mass-loss rate as a free parameter These considerations lead to 
the definition of the adopted grid of atmospheric models, sum- 
marised in Table [1] For the dusty models, we considered four 
different values for the condensation temperature of the dust 
grains (r^ = 800, 900, 1000, and 1200 K). In conclusion, for 
both groups we ran a set of 24 models; for the "Red" group we 
have 13 C-rich (two with four different Tc plus five dust-free) 



and 1 1 O-rich (two with four different Tc plus three dust-free) 
models; for the "Blue" group we have 12 (three with four differ- 
ent Tc) C-rich and 12 O-rich models. F or all stars we adopted a 
dista nce of 50 kpc (e g-, Schaefed20(j8|). an extinc t ion A u - 0.25 
mag (ISchlegel et al]|1998l) . and the lCardelli et al.l (119891) extinc- 
tion law. 
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Table 1: The model grid. For each group ("Blue" and "Red") 
and chemical composition, we list the temperatures adopted for 
models with the mass-loss rate set at zero ("0") or as a free pa- 
rameter ("F"). 



T-eff [K] 


Blue 


Red 




C-rich 


O-rich 


C-rich 


O-rich 


2600 


F 




F 




2800 






F 


F 


3000 


F 




F 




3200 





F 




F 


3600 





F 




F 


4000 













The model that yields the minimum was taken as best- 
fitting solution. Figure |2] shows - as an example - the data 
and the best-fitting models for three stars. The complete list 
of figures with the SEDs and the best-fitting models for all 
stars are available as online material, while the values are 
listed in Table 6. The parameters of best-fit models are pre- 
sented in Table |2] All stars with a Xr higher than 300 have 
been flagged as B ("Bad fits"). Most (e.g., the reddest star in 
Fig. |2]i are extremely red stars, i.e. high-amplitude long-period 
variables. The high value can then be understood consider- 
ing that our data are obtained from observations carried out at 
different epochs. Seven sources however have SEDs incompati- 
ble with AGB stars and they were therefore excluded from our 
analysis. Four (IDs 73, 130, and 133 and 165 in Tablej2j_are 
classified as candidate young stellar objects (YSO) t Woods et al.l 
1201 U iGruendl & Chull2009l). one (ID 4 3) is a Seyfert galaxy at 
z = 0.064 (IVeron-Cet tv & V6rorill2010h . one (ID 14) is a post- 
AGB star (jwoods et alJ l201lh . The last object (ID 7 5), is the 
source J050343.02-664456.7 in lGruendl &"Chul (|2009|) which is 
generically classified as an Extremely Red Object (ERO). This 
is one of the two objects not detected by VMC which were in- 
cluded in our sample using the mid-IR selection criteria. The 
classification of this source is doubtful. Gruendl et al. (20081) ar- 
gued it is an AGB star, but in the Appendix we show that its near- 
IR emission is unlikely to match the SED expected for an AGB 
star. We suggest that it may be a transition object, which left the 
AGB phase ~ 100 years ago, and entered post- AGB phase. Thus 
we excluded this star from our sample. The final sample of ob- 
jects with SED compatible with AGB stars is therefore 367. The 
photometry for all of them is given in Table 6, available as online 
material. 

The photospheric model of the best-fitting solution gives also 
a useful indication for classification of the stellar atmosphere as 
C- or O-rich. To estimate the confidence level of this classifi- 
cation, we calculated the relative difference between the two;\^^ 
values of the two best-fitting solutions obtained considering O- 
and C-rich models separately. A small difference indicates that 
the solutions obtained with C- and O-rich models are equally ac- 
ceptable; if the difference is large, the solution with the lowest 
X^ gives a noticeably better description of the data, and in this 
case the classification has a higher confidence level. After some 
tests, a value of dx^ = Ax^/Xi - 1 was adopted as discriminant. 
Stars with 5x^ < 1 were flagged as uncertain. Two tests as to 
the reliability of our classification - based on spectroscopic data 
from the literature- are provided in Sect. 14.11 

To check for any systematic effect in our modelling. Figs. 
[3] and ID show the differences between the observed magnitudes 
and those calculated from the best-fitting model in each photo- 



Table 2; Parameters of best-fitting models for all sources in our 
database. Only the first entries are shown here as an example. 
The full table is available at the CDS. The fourth column is 
the optical depth at 0.55 fim. For stars without a dusty enve- 
lope, the logarithm of the mass-loss rate (5th column), was set to 
-99.999. The 7th column indicates the C- and O-rich classifica- 
tion. Sources flagged "B" indicate stars with bad fits while "U" 
stands for uncertain classification and "S" for secure classifica- 
tion. 



id 


RA (J2000) 


DEC (J2000) 


To.55 


logM/Mo L/Lq 


class 


flag 


1 


05:02:00.86 


-65:23:12.3 


0.000 


-99.999 


3263 


O 


U 


2 


05:02:01.65 


-66:16:00.2 


0.000 


-99.999 


4847 


O 


U 


3 


05:02:02.88 


-66:46:51.6 


0.137 


-7.462 


5645 


c 


s 


4 


05:02:07.10 


-66:56:03.1 


0.019 


-8.596 


5402 


o 


s 


5 


05:02:09.01 


-66:18:57.2 


0.340 


-7.069 


5544 


c 


s 



5 



o 




0.5 
log X [/u.m] 

Fig. 2: Photometric data points (open circles) and the fitted SEDs 
(lines) for three AGB stars: a blue O-rich star (J — K^- 1 .2 mag, 
id 62), and two C-rich stars (J - Ks - 1.5 mag, id 202; and 
J-K,~4 mag, id 276). 



metric band. The variations of the distributions' spreads can be 
explained considering that AGB stars are variable sources, with 
amplitudes that are higher for redder stars and in the bluer optical 
bands. At the bluest wavelengths {B and V), the variation ampli- 
tude is maximum, as is the photometric error, since all stars - 
especially those in the red sample - are faint (many red stars are 
not detected in the B band). On the other side of the spectrum, 
due to the limited sensitivity of the Spitzer Space Telescope in 
the reddest bands, most of the bluest stars are at the detection 
limit (only 45 out of the total of 274 blue stars are detected in the 
MIPS 24 fim band). The distribution of the residuals indicates a 
good quality of the fit for the stars in the blue sample, consid- 
ering the photometric errors and the intrinsic variability of AGB 
stars. The residual distributions are in general broader for stars 
in the red sample; this is not surprising, since these are the stars 
with the highest amplitude variability, where the effects of the 
circumstellar dust are more significant. The scatter in the near- 
IR regime is similar in the red and in the blue samples, showing 
that the dust model - which mostly affects the SEDs of the red 
stars - is satisfactory. Figure |4] shows that there are systematic 
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Table 3: Results of the comparison with the classification 
by Groenewegen et al. (2009) and Groenewegen (in prep.). 
Numbers in parentheses refer only to stars with robust classi- 
fication (see text). 



-0.5 0.5 -0.5 0.5 

Fig. 4: Same as Fig. [3] but for stars with J - > 1 .5 mag. 

offsets in the reddest magnitudes. Models are always too faint in 
the IRAC 8 yum and too bright in the MIPS 24 /im bands (see 
also Fig. |2]i. This effect is negligible for stars in the blue sam- 
ple (Fig. O. The geometry of the dust shell, as well as the size 
and chemical composition of dust grains, influence the details 
of the dust emission, and different combinations of these param- 
eters could in principle produce SEDs with different slopes at 
the red end. The systematic off'sets in our models are however 
no t completely surpris ing, since a similar effect can be seen also 
in iGroenewegen et all (2009). They used a similar SED fitting 
technique, but with the addition of Spitzer spectra, giving fun- 
damental information on the dust emission features. They could 
therefore obtain much more information about the dust proper- 
ties. In this paper we assumed the dust properties that provided 
the best description of the dust emission. 

4. Results 

4. 1 . Reliability of the classification 

To test the reliability of our C-/0-rich classification, we con- 
ducted two different tests. First, we considered the cata- 
logue of C-stars in the LMC spectroscopically classified by 
iKontizas et al.l (l200lh . We found 87 stars in common with their 
sample and all are relatively blue {J - Kg < 2 mag). Our proce- 
dure correctly classified as C-stars 87% of these. If we exclude 
the 20 stars with uncertain classification, this fraction increases 
to 96%. All 54 stars with (/ - Ks > 1.5 mag) were correctly 
classified, while the success rate for the 33 bluer {J - Ks < 1.5 
mag) stars is 67%, or 84% excluding the 14 stars with uncertain 
classification. 

Second, we consider ed the stars in the catalogue of 
IGroenewegen et aP (l2009l) . complemented with other LMC and 
SMC AGB stars with Spitzer IRS mid-IR spectra (Groenewegen, 
in prepQ). The catalogue consists in 260 stars, and none of them 





Blue 




Red 








C 


O 




C 




o 




total 


46 


72 




113 




29 






(38) 


(44) 




(93) 




(22) 




correct 


93% 


72% 




97% 




76% 






(95%) 


(82%) 


(100%) 




(73%) 




Table 4: Classification of AGB stars in the VMC 8. 


.3 field. 






Blue 


Red 


All 






C 


O 


C 


O 


C 





Total 




113 


161 


93 





206 


161 


certain classification 


54 


93 


78 





132 


93 


uncertain classification 


58 


67 


10 





68 


67 


bad fit 




1 


1 


5 





6 


1 



is located in VMC field 8_3. The photometric data were there- 
fore taken from Groenewegen et al. For homogeneity with our 
analysis, only data in optical, 2MASS and Spitzer photometric 
bands were considered. We applied to the final catalogue of 260 
stars our fitting procedure and compared our classification with 
that derived by Groenewegen et al. (2009, in prep.). The results 
are summarised in Table|3] The reliability of our classification is 
confirmed to be extremely high for C-stars. In this case we ob- 
tained a 95% confidence level for classification of Blue C-rich 
stars, which is slightly hi gher than the value o btained from the 
previous test based on the Kontizas et al.l (1200 ih results, but fully 
consistent with it considering the low number of stars in the sam- 
ples used for the two tests. For O-rich stars the confidence level 
is around 75%, similar as for red stars. 

As an additional t est, we considered the seven objects in the 
SAGE-Spec sample ( IWoods et al.l 1201 ll) located in the region 
used for our study. Although this number is extremely small, the 
spectroscopic SAGE-Spec classification is in perfect agreement 
with our results. One of them is one of the YSOs and another 
one is the post-AGB star already discussed. Three sources are 
red supergiants rightly excluded by the colour selection. Finally 
two so urces are classified as C-AGB stars both by i Woods et al] 
(l20Ilh and by us. 



4.2. C- and O- rich stellar populations 

The results of our classification of the 367 AGB stars in VMC 
field 8_3 are summarised in Table |4] and in the two upper pan- 
els of Fig. [T] The location of C-and O-rich stars in the CMD 
agr ees with prediction s of AGB stellar models (e.g., see Fig. 7 
in iMarigo et al.ll2008l) : O-rich stars with no -or small- dusty 
envelopes are located on a nearly vertical blue plu me, while C- 
rich stars populate the reddest part of the CMD. In lMarigo et al.l 



' Compared to IGroenewegen et al] ( j2009h . AGB and RSG stars 
from the following IRS programs were added: 30788 (PI: Sahai , 



no dedicated publication, 3 spectra are discussed in iBuchanan et all 
20 0l); 40159 (PI: Kemper, the SAGE-SPEC program described i n 
,Kemr)eretalJ 120101) : 40650 (PI: Looney, see iGruendl et id] l2008h : 
50167 (PI: Cl ayton, no publicat ion); 50240 (PI: Sloan, the SMC-SPEC 
program, see ISloan et al.ll2010D . 
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Fig. 3: Histograms of magnitude residuals for the data in all wave-bands, for O-rich {left panel) and C-rich {right panel) stars with 
J - Ks < 1.5 mag. The histograms are normalised to the total number of stars detected in each wave-band, labelled in each plot. 



(120081) this is ascribed to the cool Tefi caused by changes in 
molecular opacities as the third dredge-up events increase the 
photospheric C/O ratio. 

The most remarkable fact is that all red stars are found to 
be C-stars. We cannot exclude, in principle, that some of them 
could be misclassified O-rich. The results of our tests presented 
in Sect. 14.11 showed that the probability of an incorrect classi- 
fication is around 25%. The presence of more than one O-rich 
dust-enshrouded star can therefore be excluded at a ~ 94% con- 
fidence level. 

A very small number - if any - of red O-rich stars was 
expected a priori since dust-enshrouded O-rich stars are rare: 
low mass O-rich AGB stars do not reach high-mass rate ( > 
10^^ Mq yr '). On the other hand, the most massive O-rich AGB 
stars are intrinsically rare objects, because of the shape of the 
stellar initial mass function, which is dominated by low-mass 
stars in a power-law fashion. Dusty O-rich stars are therefore 
expected to be found preferentially in regions with a high star- 
formation rate at recent epochs (< 1 Gyr). The VMC 8_3 field, 
used for this preliminary study, is located in the outer regions of 
the L MC, where the rece nt star-formation rate is not particularly 
high dRubele et al.ll2011h . 

All the brightest stars in the blue verti cal plume at Ks < 1 
mag are classified as O-rich, as expected iMarigo et alJ (l2008h . 
That region of the CMD is in fact populated by the massive O- 
rich AGB stars in the hot bottom burning phase. However note 
that we cannot exclude that some of them could in principle be 
core He-burning red supergiants. 

4.3. Bolometric magnitudes 

The bolometric correction BCk - M\,o\ - Mk, obtained from the 
total luminosity and the /T^-band magnitude (corrected for ex- 
tinction and distance modulus) is shown in Fig.|5]as a function 
of J - Ks colour. To extend our data point s to redder colours, 
we included also LMC and SMC data from lGroenewegen et alJ 



o 
m 




Fig. 5: Bolometric correction as a function of near-IR colour O- 
rich stars are shown as blue open squares and C-rich stars as 
filled red c ircles. Crosses are data for C-stars in the LMC and 
SMC from iGroenewegen et alj ( l2009h . T he open black circle s 
show the relation for Gal actic C-stars from Bergeat e t al.l (l2002h . 
the lines are taken from iKerschbaum et al.i (.2010.) : the curved 
line is for C-stars, while the two straight lines define the region 
of O-rich stars. 



20091). Our data are fully consistent with ot her literature data 



Bergeat et alJl2002l : IKerschbaum et al.ll2010l) . The only (possi 
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en 

O 







O 




Fig. 6: Bolometric luminosity function of AGB stars. The filled 
histogram and the dashed line show the LF obtained from C- 
rich and O-rich stars, respectively. The inner panel shows the LF 
for dust-free (dotted) and dusty (solid) AGB C-stars. The solid 
curve is a Gaussian fit to the histogram of dusty stars. 



Fig. 7: Mass-loss rates as a function of IR colours. Symbols are 
the same as in Fig.|5] The grey open histogram shows the colour 
distribution of dust-free AGB stars; the red histogram is the dis- 
tribution of C-stars alone. The solid line is th e best fit to all C- 
stars, including the lGroenewegen et al.l ( 2009|) data. The dash ed 
line in the upper panel is the relation of iMatsuura et al.l (l2009h . 



ble) exception could be represented by the extremely red C-stars 
at / - Ks > 2.5 mag. For these st ars the BCs obt ained from both 
our data and t hose of [Or oenew egen et al.l (l2009h are located 0.3 
mag below the lKerschbaum et alT ( 2010h relation. Note however 



that this r elation is not well defined a t J-Ks > 2.5 mag, since the 
sample of iKerschbaum et all (l201d) contains only few red stars, 
which exhibit a significant scatter and lie - on average - below 
the best-fitting relation. We can therefore conclude that the bolo- 
metric magnitudes obtained from the SED fits are in agreement 
with literature data at a level of - at worst - few tenths of mag- 
nitude. 

The bolometric LF function obtained for all AGB stars in 
our sample is shown in Fig. |6] The cut at the faint end, around 
Mhoi = -3.7 mag, is the consequence of the selection cut for 
stars fainter than the TRGB. The overall LF peaks at Mboi - -4.5 
mag and drops at magnitudes brighter than Mboi ~ -5.0 mag. 
There are very few stars brighter than Mboi = -6. and none 
above the classical limit for AGB stars, Mboi - -7.1 dPaczvnskil 
Il97lh . corresponding to a Chandrasekhar-limit core-mass be- 
fore core He ignition. The most massive AGB stars, under- 
going "hot bottom burning", may h ave a luminosity brighter 
than the classical AGB limit (s ee, e.g.. lGroenewegen et al.l2009t 
iGarcfa-Hernandez et al.ll2009l) . but their number is expected to 
be extremely small. 

The LF of C-rich stars shows a bimodality, which can be 
understood by considering separately the LFs for stars whose 
best-fitting solution is dust-free and and stars with dusty en- 
velopes. Dust-free stars show a decreasing LF, truncated by 
the selection c ut. Dusty C-sta r s are typically brighter - as al- 
ready noted by Ivan Loon et al ] (119991) - and their LF is well de- 
scribed by a Gaussian function centred at Mboi - -4.84 mag and 
with cr - 0.40 mag. The peak of the LF is in perfect agree- 
ment with other observations of LMC C-stars and fully con- 



sistent with the expectations of theoretical models (see, e.g., 
iGroenewegen & de Jondfiggl . Some of the faintest AGB stars 
classified as C-rich could actually be misclassified O-rich stars, 
since our classification procedure is less precise for faint and 
blue stars, but in any case, our result seems to indicate the pres- 
ence of a non negligible population of AGB stars fainter than 
the TRGB characterised by very low - if any - mass-loss rates. 
This is in agreement with the prediction of recent evolutionary 
models, in contra st with the paucity of faint C stars in spec- 
troscopic surveys ([Kontizas et al.ll2001l : lGulUeuszik et al.ll2008l : 



iMarigo et al.ll2008l) . We will reserve a more detailed discussion 
of this issue for a future paper, presenting a comparison of our 
data with the predictions of theoretical evolutionary models. 

Recently, ISrinivasan et al.l (1201 ll) presented a C-star LF, re- 
vising their earUer work ( Srinivasan et al.l l2009h . The general 
shape and the position of the peak in the revised LF is in good 
agreement with our result. 

4.4. Mass-loss rates 

Figure |7] shows mass-loss rates as a function of IR colours. The 
data points lie on a well defined sequence, fully compatible with 
data published in Groenewegen et al. ( 2009), showing that our 
fitting technique, based only on photometric observations, pro- 
vides reliable mass-loss rates. Also in this case, the main limi- 
tation of our results is due to the poor statistics for the reddest 
dusty stars. We have no O-rich stars with noticeable mass-loss. 
The number of dusty C-stars is still very small, but the sequence 
of the few data points redder than J - Ks ~ 2 mag exhibits a 
remarkably low dispersion. From a least-squares fit to all data. 
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including those of iGroenewegen et al.l (l2009l) . we obtained the 
following relation for C-stars: 



logM : 



-15.42 



■2.70 



(7 -/*:,) + 2.10 
for l.O < J - Ks < 5.0 mag; and 



logM : 



-12.78 



([3.6] - [8.0]) -t- 2.49 



2.63 



(2) 



(3) 



for 0.0 < [3.6] - [8.0] < 4.5 mag. For the bluest stars in our 
sample the mass-loss rates are extremely low, at the detection 
limit of our method. At J - Ks ~ 1 .0 mag, the best-fitting model 
was dust-free for most stars; for the others, the mass-loss rates 
are less than ~ lO^^Moyr"'. A conservative approach would 
suggest that the measured values for the bluest stars must there- 
fore be considered upper limits. The blue applicability limits of 
Eqs. (|2]i and O) were based on this consideration and on the 
colour distributions for dust-free stars in Fig. |7] On the other 
hand, the definition of the red limit is somewhat arbitrary and 
it was set by the colour of the reddest stars in our sample. 
The up per panel of Fig. Ill co mpares our results with the rela- 
tion of iM atsuura et al.' ("2009^, obtained using all measures of 
IGroenewegen et al., (.2009i) . The difference for stars bluer than 
[3.6]-[8.0]= 2 mag is main ly due to the small numbe r of stars 
in this colour range in the Groenewegen et al.l (|2009|) sample. 
The relation of IMatsuura et al. ( 20091) is therefore weakly con- 
strained for the bluest stars. In addition, their relation is derived 
using the colours obtained from the best-fitting SEDs, rather than 
the observed counterparts, as in our case. The difference is there- 
fore related to the systematic offsets already discussed in Sect. 
[3] We conclude that our relation is more reliable when mass-loss 
rates are estimated from observed photometry. 



5. Discussion and Conclusion 

In this paper we used the first data from the ESO public survey 
VMC to obtain optical to mid-IR SEDs for a sample of 367 AGB 
stars candidates. Our selection criteria were defined to minimise 
the selection bias. Sources with colours as blue as the RGB and 
fainter than Kg = 12.0 mag -which roughly corresponds to the 
tip of the RGB in the LMC- were excluded since this region of 
the CMD is dominated by the RGB and it is nearly impossible to 
disentangle RGB and AGB stellar populations from photometric 
data. 

Luminosities and mass-loss rates were obtained by fitting the 
SEDs with a dust radiative transfer model. Mass-loss rates were 
obtained assuming a constant expansion velocity and dust-to-gas 
ratio with values similar to Galactic AGB stars. 

The cumulative distribution of mass-loss rates for O- and 
C-rich stars is shown in Fig. [8] The total mass-loss rate in 
the 1.42 deg- area encompassed by our analysis is 6.4 x 
lO^^Moyr"'. This value is lower than that expected from the 
average mass-loss rate for the external regions of the LMC 
found bv iMatsuura et all (l2009l) . which is 1.1 x lO^'^Moyr ' - 
scaled to a 1 .42 deg^ area. We note however that the ERO ob- 
ject discussed in the A ppendix was excluded from our analysis. 
IMatsuura et al.l (l2009l) considered it to be an AGB and estimated 
for it M = 6.4 X 10"^ Mq yr This object alone could explain 
the difference between the integrated mass-loss rates. 

Table |5] shows the star counts and the integrated mass-loss 
rates for stars in different mass-loss ranges, toget her with the star 
counts for the external region of the LMC from IMatsuura et al.l 




Fig. 8: Cumulative distribution of the mass-loss rates for O- 
{open squares) and C-rich (filled dots) AGB stars. The horizon- 
tal line shows the averag e mass-loss rates of the external regions 
of the LMC obtained by IMatsuura et al.1 (l2009l) . rescaled to the 
area encompassed by our analysis. 

Table 5: Distribution of mass-losing O- and C-rich sta r s. The 
second column shows the star counts of IMatsuura et al.l (l2009l) 
for C-stars in the external region of the LMC rescaled to the area 
encompassed by our analysis. The last two columns are the star 
counts and total mass-loss rates obtained from this study. 



C-rich 





M range 


Nm09 


N 


Mtot 




Me yr' 






10"^ Mo yr"' 




< 1 X 10-^ 


9.1 


102 


1.5 


1 X 10 


* < M < 3 x 10"^ 


8.1 


6 


1.1 


3x 10 


<M <6x 10-" 


4.4 


3 


1.3 


6x 10 


" <M <lx 10-^ 


1.8 








1 X 10 


<M <3x 10-^ 


1.8 


2 


2.5 


3x 10 


<M <6x IQ-^ 


0.5 










> 6 X 10--^ 


0.2 










Total: 






6.4 


O-rich 




M range 


Nm09 


A' 


M-poT 




< 1 X 10-^ 




65 


0.05 




Total: 






0.05 



( 120091) . Since lMatsuura et alj ( |2009|) considered only C-stars, O- 
rich stars are not taken into account in Table [5] The star counts 
are in good agreement, and the only significant discrepancy is re- 
lated to stars with lo w mass-loss rate s . This is due to the selection 
criteria adopted by IMatsuura et al.l (l2009l) . biased towards red 
stars. The 102 stars with M < 10"^ Mq yr"' contribute to ~ 24% 
of the total mass-loss rate, and the two stars with higher mass- 
loss rates contribute up to 38 % of the total. The fact that the inte- 
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grated mass-loss rate is strongly dominated by the most extreme 
stars, prevents us from drawing any conclusion on the overall 
mass-loss rate in the LMC from the relatively small number of 
stars in our sample. The highest values of the ma ss-loss rates in 
the L MC can be as high as ~ 10 "^Moyr ' (e.g. jGruendl et al.l 
l2008h . To reach a more robust estimate we therefore need to ex- 
tend our analysis to a much larger area. The 65 mass-losing O- 
rich stars have an integrated mass-loss rate of 4.5 x 10"^ Mq yr"' , 
which is negligible compared to the C-stars. Nevertheless, this is 
not true in general, since extreme O-rich AGB stars can have 
mass-loss rates as hi^h as the integrated value of our sample 
jvan Loon et a l. 2005V Groenewegen et al.ll2009l) . 

The main limitations of the results reported in this paper are 
related to the relatively small number of dust-enshrouded AGB 
stars in our sample. The main aim of our work was to explore 
the confidence level of the mass-loss rates and bolometric mag- 
nitude measures derived from photometric SEDs constructed by 
combining VMC data with optical MCPS, near-lR 2MASS and 
mid-IR Spitzer data, with a view to applying our method to the 
complete sample of AGB stars that will be observed by VMC 
in the Magellanic system. We demonstrated the reliability of the 
measures and classification of the O- and C-rich AGB stellar 
populations. A more complete picture about the mass-loss re- 
turn and the process of enrichment of the ISM will be obtained 
from forthcoming VMC observations. VMC is scheduled to reg- 
ularly carry out observations which, at the end of the 5 year sur- 
vey, will provide photometry for two orders of magnitude more 
AGB stars. Our results are nevertheless important, because our 
database already contains a sufficiently large number of stars 
with colours up to / - Ks - 2 mag, which represent the bulk 
of the AGB population. 
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Appendix A: J050343.02-664456.7 

This section is dedicated to a more detailed analysis of the ex- 
tremely red source J050343. 02-664456. 7. It is one of the red- 
dest sources in our sample, and one of two objects detected only 
in the mid-lR, with no counterpart in the VMC catalogue - the 
other one being a YSO. 

J050343 .02-66445 6.7 is one of the 13 sources identified by 
iGruendl & Chul (l2009h and generically classified -on the basis 
of their photometric SED- as EROs. Subsequent Spitzer/IRS 
spectroscopic observation s of seven of them r evealed that they 
are extreme carbon stars (iGruendl et al.ll2068h . The conclusion 
is based on the detection of SiC and C2H2 absorption and on the 
presence of a broad MgS feature in two cases. 

Assuming that these objects are AGB stars, the derived mass- 
loss rates are hi gher than those of an y known carbon-rich AGB 
star in the LMC. IGruendl et akl (l2008l) estimated a mass-loss rate 
M^l.Sx IQ-'^ Mo yr-i for J050343.02-664 456.7. From empir- 
ical relations based on mid-lR colours, iMatsuura et al ] (I2009h 
found a lower value, viz. 6.4 x 10"^ Moyr This object alone 
provides a contribution equal -or bigger- than the integrated 
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Fig.A.l: 30" X 30" J (upper panel) and (lower panel) 
VMC images centred on the Extremely Red Object J050343.02- 
664456. 



mass-loss rate of all stars in our sample. In the following we 
present some indications that J050343. 02-664456. 7 may not be 
a normal AGB star 

J050343. 02-664456.7 is not present in the VMC catalogue 
because it is detected only in the /(Tj-band. It is in fact completely 
invisible in the / band as well as in the bluer Y band, as shown 
in Fig. I A. II The vmcSource ca talogue is built wi th sources de- 
tected in the three VMC bands dCioni et al. l l20Tlh and therefore 
J050343.02-664456.7 is not included. The K, magnitude was 
obtained from the vmcDetection table, i.e. the catalogue cor- 
responding to individual observations. Its value is Ks = 18.75 
mag. 
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Fig. A. 2: Photometric data points and standard AGB SED model 
for the Extremely Red Objects J050343. 02-664456 (diamonds 
and solid line). In the upper panel they are compared with data 
and SED standard model for J055026.08-695603. 1 ((circles and 
dotted line). The lower panel show the post-AGB model with a 
detached dusty shell (dotted line). The arrows show the photo- 
metric detection limits. 



The photometric data points and the best-fitting model of 
J050343. 02-664456. 7 obtained using the standard procedure de- 
scribed in Sect.|3]are presented in Fig. lA.2l The model substan- 
tially underestimates the flux in the bluest part of the SED, in the 
Ks and in the IRAC 3.6 yum bands. The predicted Ks magnitude 
is ~ 9 mag fainter than the observe d one. We tested our m odel 
also on the other EROs identified bv lGruendl & Chul (l2009l) -aU 
located outside the 8_3 VMC tile and hence not included in our 
sample- and found that in most cases our AGB models could de- 
scribe quite well the observed SED. As an example, in Fig. lA.2l 
we show the photometric data and our best-fitting SED model for 
J055026.08695603.1, an EROs with a SED similar to the one of 
J050343. 02-664456. 7. In this case flie model SED seems to gen- 
erally better reproduce the observed data. We are therefore led 
to consider the possibiUty that J050343. 02-664456. 7 could be 
something different from an AGB star. 

The SED of J050343. 02-664456.7 could be explained as- 
suming that this is a post-AGB star with a detached shell. In 
this case the flux in the near-IR would be due to the contribution 
of the central star. In our case this contribution is extremely low, 
and produces only a flatter SED in th e near-IR, rather than a sec- 
ondary peak in the near-IR (see, e.g.. lLagadec et al ] l2011h . This 
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implies that the central star stopped l osing mass -i.e. it ended t he 
AGB phase- very recently (see, e.g.. Ivan der Veen et al ] ll989h . 

To explore the post-AGB hypothesis, we used our models 
setting the condensation temperature of the dust grains as a free 
parameter The resulting best-fitting model is shown in the lower 
panel of Fig. IA.2l as a dotted line. It is compatible with the up- 
per limit for optical and near-lR magnitudes corresponding to 
the MCPS and VMC detection limit and it shows a much bet- 
ter agreement with the observed SED than the standard AGB 
model. In the best-fitting model the condensation temperature 
is Tc - 330 K, corresponding to an inner radius of the dusty 
shell of ~ 5 X IQ'^Rq or ~ 3.5 x 10'" km. This value is ex- 
tremely small, just 10 times bigger than the inner radius of the 
dusty shell obtained for the standard model. This implies that 
the mass-loss production in J050343. 02-664456. 7 has dropped 
to zero extremely recently -of the order of ~ 100 years ago as- 
suming a shell expansion velocity of 10 km s"' . 

To conclude, we showed that the observed SED of 
J050343.02-664456.7 is not fully compatible with a standard 
model for AGB stars, showing a flux excess at wavelengths 
shorter than ~ 4 /^m. We proposed that this source could be- 
long to the -rare- class of objects in transition from the AGB to 
the planetary nebula stage. Further evidence supporting this idea 
could be fou nd in the fact that some of the SEDs of the 13 EROs 
identified bv lGruendl&Chul (l2009h show hints of a secondary 
peak at near-IR wavelengths. We hence point out that it may be 
not so straightforward to classify all EROs as AGB carbon stars. 



